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We consider here scalar and electromagnetic perturbations for the Vaidya metric in double-null 
coordinates. Such an approach allows one to go a step further in the analysis of quasinormal 
modes for time-dependent spacetimes. Some recent results are refined, and a new non-stationary 
behavior corresponding to some sort of inertia for quasinormal modes is identified. Our conclusions 
can enlighten some aspects of the wave scattering by black holes undergoing some mass accretion 
processes. 
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I. INTRODUCTION 

The quasinormal modes analysis is a paradigm for 
the study of gravitational excitations of black holes. 
For a comprehensive review on the vast previous lit- 
erature, see P, Q- We stress that the quasinormal 
modes (QNM) of black holes are also experimentally 
relevant since they could, in principle, be detected in 
gravitational waves detector such as LISA (see, for in- 
stance, |3| and references therein). The robustness of 
QNM analysis has been established by considering sev- 
eral kinds of generalizations, as, for instance, relax- 
ing the asymptotic flatness condition^, JM Q or con- 
sidering time-dependent situations 0, M> 01 • The case 
of asymptotic Anti-de-Sitter black holes has deserved 
special attention due to its relevance to the AdS/CFT 
conj e c t ur e [13 111 111 111 111 111 111 [11 111 lii S 111 

The analysis of QNM for time-dependent situations is 
of special physical interest since it is expect that a black 
hole can enlarge its mass by some accretion process, or 
even lose mass by some other process, including Hawking 
radiation. The late time tail for a Klein Gordon field 
under a time-dependent potential was first considered in 
13, where some influence of the temporal dependence of 
the potential over the characteristic decaying tails was 
reported. Such a kind of time-dependent potential arises 
naturally when the Vaidya metric is considered .8, 9] . The 
Vaidya metric, which in radiation coordinates (w, r, 6, (j)) 
has the form 



2m{w) 



2cdrdw + r'^dn'^, (1) 



where dQ'^ = dO'^+sm'^ 



c = ±1, is a solution of Ein- 



stein's equations with spherical symmetry in the eikonal 
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approximation to a radial flow of unpolarized radiation. 
For the case of an ingoing radial flow, c = 1 and m(w) 
is a monotone increasing mass function in the advanced 
time w, while c = — 1 corresponds to an outgoing radial 
flow, with m{w) being in this case a monotone decreasing 
mass function in the retarded time w. The metric ^ is 
the starting point for the QNM analysis of varying mass 
black-holes. 

In 9], massless scalar fields were studied on an elec- 
trically charged version of the Vaidya metric. Basically, 
two kinds of continuously varying mass functions were 
considered: 



mi w <wi, 
m{w) = \ mi/(w) wi <w < W2, 

1712 W >W2, 



(2) 



with f{w) = (1 — Xw) and f{w) = exp(— au;), called, 
respectively, linear and exponential models. General- 
ized tortoise coordinates were introduced, and a stan- 
dard numerical analysis was done. The conclusion was 
that, as a first approximation, the QNM for a station- 
ary Reissner-Nordstrom black hole with mass m{w) and 
charge qi^w) = qom(w) are still valid, as if some sta- 
tionary adiabatic regime was indeed governing the QNM 
dynamics for time-dependent spacetimes. In principle, 
one should not expect such stationary behavior for very 
rapid varying mass functions, however the analysis of [3 
was not able to identify any breaking of such an adiabatic 
regime. We notice also that both linear and exponential 
models, inspired by some known exact solutions^ for 
which generalized tortoise coordinates could be explic- 
itly constructed, would hardly correspond to physically 
realistic cases. Both cases have C°-class mass functions, 
implying the existence of some infinitesimal shell distri- 
butions of matters for w = wq and w — wi whose inter- 
pretation and role are still unclear. 

We attack here the QNM problem by considering the 
Vaidya metric in double- null coordinates 24] . The main 
advantage of such an approach is the possibility of con- 
sidering any (monotone) mass function, allowing, for in- 
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stance, the analysis of smooth mass functions that could 
correspond to physically more relevant situations, free 
of obscure infinitesimal matter shells. Furthermore, the 
use of double- null coordinates has improved considerably 
the overall precision of the numerical analysis, allowing 
us to identify the breakdown of the adiabatic regime, 
characterized by the appearing of a non-stationary in- 
ertial effect for the QNM in the case of rapid varying 
mass functions. Once the stationary regime is reached, 
the standard QNM results hold, reinforcing, once more, 
the robustness of the QNM analysis. Our results can 
be used as a first approximation to describe the propa- 
gation of small perturbations around astrophysically re- 
alistic situations where black holes undergo some mass 
accretion processes. We notice that the non-stationary 
inertial effect presented here is in perfect agreement with 
the non-linear analysis preformed in pH , which reported 
extra rcdshift effects on the QNM frequencies due to the 
growth of the effective black hole mass. 

We organize this paper as follows. The next section 
is devoted to a brief introduction to the semi-analytical 
approach|24j for the Vaidya metric in double-null coordi- 
nates. Section III presents the main issues and results of 
our numerical analysis. The last section is left to some 
closing remarks. 



II. THE VAIDYA METRIC IN DOUBLE-NULL 
COORDINATES 

Double-null coordinates are specially suitable for time- 
evolution problems as the QNM analysis. However, the 
difhcultics of constructing double-null coordinates for 
non-stationary spacetime are well known. The Vaidya 
metric is a typical example (see p^ ] and references 
therein). Indeed, it is known that the problem of con- 
structing double-null coordinates for generic mass func- 
tions is not analytical soluble in general _2_3:j. The semi- 
analytical approach proposed in |2J| is the starting point 
for our analysis here. It consists, basically, in consider- 
ing the Vaidya metric in double-null coordinates ab initio, 
avoiding the need of constructing any coordinate trans- 
formation. The spherically symmetric line element in 
double-null coordinates is 



-2f{u, v)du dv + {u, v)dn'^ , 



(3) 



where f{u, v) and r(u, v) are smooth non vanishing func- 
tions. The energy-momentum tensor of a unidirectional 
radial flow of unpolarized radiation in the eikonal approx- 
imation is given by 



Tab = ^h{u,v)kakb, 



(4) 



where ka is a radial null vector. The Einstein's equa- 
tions for the case of a flow in the u-direction can then be 
reduced to the following set of equations ;2S.i24| 



dyr{u,v) 
h{u, v) 



-B{v) 1 



2to(u) 



r(u, v) 



^B{v)m' (v) 

arbitrar; 
;y condil 

B{v)m'{v) < 0, 



(6) 
(7) 



where B{v) and m(w) are arbitrary functions obeying, 
according to the weak energy condition, 

(8) 

where the prime denotes the derivative with respect to 
V. The solution of ©-0 will correspond to the Vaidya 
metric in double- null coordinates, as one can interpret 
from (O and l@J. For m'{v) ^ 0, the choice 



B 



--sign (to') 



(9) 



allows one to interpret m(v) as the mass of the solution 
and V as the proper time as measured in the rest frame 
at infinity for the asymptotically flat case [23, 13- Note 
that if the weak energy condition (jSJ holds, the function 
m{v) is monotone, implying that the radial flow must 
be ingoing or outgoing for all v. It is not possible, for 
instance, to have "oscillating" mass functions m{v). 

The semi-analytical approach of 2^ consists in a strat- 
egy to construct numerically the functions /(u, w), r(u, v) 
and h{u, v) from the Eq. ©-Q, and to infer the under- 
lying causal structure. Eq. 10 along u constant is a 
first order ordinary differential equation in v. One can 
evaluate the function r{u,v) in any point by solving the 
u-initial value problem knowing r{u,vo). The trivial ex- 
ample of Minkowski spacetime (to = 0), for instance, can 
be obtained .24| by choosing r{u, 0) = u/2. Once we have 
r(u,v), we can evaluate f{u,v) and h{u,v) from Q and 
0. Further details of the method can be found in [23| . 
We consider here the following smooth mass function 



m{v) = TOi -|- 



7712 — TOl 



[1 + t&nh p{v - vi)] (10) 



where p and vi are also constant parameters. For sake 
of comparison with the results of 9], we also consider a 
linear model. Fig. 1 depicts the causal structure corre- 
sponding to the hyperbolic mass function (|10() obtained 
from the semi-analytical approach. For positive p and 
TO2 > TOi, for instance, it represents a black hole with 
mass TOi receiving a radial flux of radiation and, con- 
sequently, enlarging continuously its mass until reaching 
TO2. The choice of the initial condition for solving © 
is a rather subtle issue^- For our purposes here, we 
note only that demanding dur{u,vo) 7^ is sufficient to 
guarantee that, for any m{v), the underlying spacetime 
causal structure does not depend on the initial condition 
r{u,vo). For constant to and B = —1/2, Eq. ® can be 
easily integrated, leading to 



r{u, v) + 2m In |r(u, v) — 2m\ 



Pin), (11) 



f{u,v) = 2B{v)dur{u,v), 



(5) 



where P{u) is an arbitrary function of u. It is also possi- 
ble to solve Eq. ijHJl analytically for the linear and expo- 
nential mass functions [23]. Eq. Hll|) is our reference to 
choose the initial conditions r(w, vq). 
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FIG. 1: Lines of constant u for the solution of Q with the 
hyperbolic mass function dU, with m2 > mi and p > 0. 
All solutions in the region below the line r = 2m{v) (the 
apparent horizon) have r' < 0, where the prime denotes the 
derivative with respect to v. Any solution that enters into this 
region will reach the singularity at r = with finite v. On 
the other hand, solutions confined to the r' > region always 
escape from the singularity and reach T"*". In the present case, 
there exists an event horizon (the dashed line on the inserted 
conformal diagram) close to the solutions rniv)- 



III. SCALAR AND ELECTROMAGNETIC 
PERTURBATIONS 



the values ip{umax,v) are extracted, where Umax is the 
maximum value of u on the numerical grid. Taking sufB- 
ciently large Umax, we have good approximations for the 
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FIG. 2: Values of the electromagnetic perturbations 
'4^{umax,v) with £ = 2 for a Schwarzschild black hole with 
mass m — 0.5 and for the Vaidya metric with the hyperbolic 
mass function IIUL where Vi = 75, p — 0.08, mi — 0.5, and 
m2 — 0.65 (Hyperbolic 3). One can clearly appreciate, for 
the time-dependent case, the slowing down of the oscillation 
frequency and damping taking place after v\. 



In the coordinate system for any desired form of 
the mass function m{v), the equations for scalar and 
electromagnetic perturbations can be easily put in the 
formQ 



+ V{u,v)f{u,v)ij = 0, 



dudv 

where the potential V(u, v) is given by 



V{u,v) = 



£{£+!) 



m(v) 



2r'^(u, v) r^(u, v) 



(12) 



(13) 



where cr = 1 and a ~ correspond, respectively, to the 
scalar and to the electromagnetic cases. For a given mass 
function m{v), one first evaluates the function f{u,v) 
and r(M, v) with the semi-analytical approach of the pre- 
vious section, and then the characteristic problem corre- 
sponding to (|12|) can be solved with the usual second or- 
der characteristic algorithm^, where the initial data are 
specified along the two null surfaces u = uq and v — vq. 
Since the basic aspects of the field decay are indepen- 
dent of the initial conditions (this fact is confirmed by 
our simulations), we use the Gaussian initial condition 



ip{u = Mo, v) = exp 



{v ~ Vcf 



2a2 



(14) 



and ip{u,v = vq) = const. Our typical numerical grid 
is large enough to assure that we can set effectively this 
last constant to zero. After the integration is completed, 



wave function at the event horizon. Figure 2 presents 
an example of ip{umax,v) data for a hyperbolic increas- 
ing mass function and its comparison with the case of 
a Schwarzschild black hole. Similar results hold if we 
extract the data for other values of u. All the analysis 
presented here corresponds to the data extracted on the 
horizon Umax- Note that the choice of Umax is a matter 
of numerical convenience. Since the QNM corresponds to 
eigenstates of an effective Schroedinger equationQ, the 
associated complex eigenvalues can be read out anywhere 
once the asymptotic regime is attained. 

From the function ip{umax,v), one can infer the char- 
acteristic frequencies lo of the damped oscillating modes. 
Since the frequencies w are themselves time-dependent 
(see Fig. 2), for a given w, uj{v) is determined locally 
by a non-linear x^-fitting by using some damped cycles 
of "0 around v. We perform an exhaustive numerical 
analysis for many different mass functions m{v). Fig- 
ures 3 and 4 present the real (wr) and imaginary (wi) 
part of the frequencies of the perturbations as function 
of V. They correspond, respectively, to some decreas- 
ing and increasing mass functions, and are typical for all 
values of £, all kinds of perturbations (scalar and elec- 
tromagnetic) and different initial conditions. The imagi- 
nary part LUi exhibits similar time-dependence, although 
its values are typically known with lower precision than 
tjR. This is not, however, due to numerical errors in 
the integration procedure. In fact, the convergence of 
the second order characteristic algorithm is quite good. 
For instance, the differences in the calculated frequen- 
cies of Fig. 3 and 4 obtained by using as integration 
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unnoticed by the QNM analysis in radiation coordinates 
performed in Q • Analogous conclusions hold for the case 
of increasing m(w) (Fig. 4). In fact, since the increas- 
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FIG. 3: The real {ujr) part of the frequency for scalar pertur- 
bations with ^ = 2 as function of v for decreasing linear and 
hyperbolic mass functions. For all the cases, mi — 0.5 and 
m2 = 0.35. Linear 1: vi — 60, «2 = 90. Linear 2: vi — 74.5, 
V2 ~ 75.5. Hyperbolic 1: iii = 75, p — 0.08. Hyperbolic 2: 
ui = 75, p = 0.8. The imaginary part cui has similar behavior 
as for the increasing mass case (Fig. 4), see the text. 



steps Au = Av = 0.2,0.1, and 0.05, are smaller than 
the size of the points used in the figures, and no differ- 
ence is detected for even smaller steps. We credit the 
irregularities observed for uij to the local x^-fitting de- 
scribed above. For subcritical damped oscillations as the 
QNM considered here, the oscillation frequencies cjr can 
be easily determined from a few damped cycles, while for 
the damping term uji one typically needs many more cy- 
cles in order to get a similar precision. By taking a large 
number of cycles we will tend to smear the calculated 
frequencies, leading to some kind of average and not to 
the desired instantaneous values. We opt, then, to take 
as few as possible cycles, paying the price of having some 
irregularities for wj. 

Let us take for a close analysis, for instance, the case 
of decreasing m{v) (Fig. 3). For the rapid varying cases 
(Linear 2 and Hyperbolic 2) , one clearly sees the inertial 
effect for near v = 75. The function lur^v) does not 
follow the track corresponding to m~^{v), as one would 
expect for a stationary adiabatic regime, and as it re- 
ally does for the Hyperbolic 1 case. After the rapid in- 
creasing phase, acts as it would have some intrinsic 
inertia, reaching a maximum value that is bigger than 
(x'r(oo), implying the relaxation corresponding to the re- 
gion with uj{^{v) < for v > 75. We also notice that, 
for the rapid varying case, one could not detect sensible 
differences between the smooth hyperbolic case and the 
C° linear one. Nevertheless, one can see that the slower 
linear case (Linear 1) exhibits some inertial effects close 
the the second matching point of m{v). In all other re- 
gions, it follows the track corresponding to m~^{v). No 
appreciable difference in the transients of scalar and elec- 
tromagnetic perturbations was detected. As we have al- 
ready mentioned, such transient inertial behavior passed 
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FIG. 4: The real {uok) and imaginary (lji) part of the fre- 
quency for electromagnetic perturbations with ^ = 2 as func- 
tion of V for increasing linear and hyperbolic mass functions. 
For all the cases, mi = 0.5 and m2 = 0.65. Linear 3: «i = 60, 
V2 ~ 90. Linear 4: vi = 74.5, «2 ~ 75.5. Hyperbolic 3: 
vi = 75, p = 0.08. Hyperbolic 4: «i = 75, p = 0.8. 

ing mass functions considered here can be obtained from 
the decreasing ones by a time-reversal m{v) m{—v), 
the underlying causal structures are equivalent^], and 
the corresponding QNM are also related. Under a time- 
reversal, the graph of (and of lui) must be reflected 
on a horizontal axis. 

We could infer from our numerical simulations the 
situation corresponding to the onset of the QNM non- 
stationary inertial behavior. The deviation from the 
stationary regime is measured by the second derivative 
of the mass function m"{v). Heuristically, one should 
expect the appearing of non-stationary behavior when 
|l/m"| is smaller than a certain characteristic time of the 
system, what should prevent the system to relax into an 
adiabatic regime. There are two characteristic times as- 
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sociated to the QNM of black holes: the oscillation period 
27r/a;R and the damping time |l/wi|. The onset of the in- 
ertial behavior is associated to the second one. We verify 
appreciable deviations from the stationary regime when- 
ever |to"(w)| is of the same order (or larger) than \lui\. 
For instance, for Hyperbolic 2 and Hyperbolic 4 data set 
(corresponding to the inertial behavior depicted in Figs. 
3 and 4) we have \m'^^^/ uj f'^'^^l w 25%, while for Hyper- 
bolic 1 and Hyperbolic 3 (the stationary behavior) such 
ratio is 100 times smaller. Estimating the magnitude of 
the inertial effect from our numerical simulations with- 
out an approximated analytical model seems to be much 
harder. Again heuristically, one expects that the magni- 
tude of the effect be proportional to |to"|t, where t is the 
time interval along which \m"\ ^ \uji\. We could check 
that the relative variations in frequencies during the in- 
ertial behavior for the smooth case are always limited by 
the ratio \m" /lui\. Further analytical work is certainly 
necessary to enlighten this point. We finish by notic- 
ing that for the C" linear case, is always large, 
increasing indeed with 1/Av, and hence non-stationary 
inertial behavior is reported in all simulations using C'^ 
linear data set. Since r is very small in such cases (of the 
same order of A?;), the magnitude of the inertial effects 
are also typically small. 

IV. FINAL REMARKS 

All situations we considered in the present work involve 
mass functions corresponding to an initial black hole with 
mass mi undergoing some accretion/losing mass pro- 
cess and ending with a mass m2. Such an "asymptotic 
Schwarzschild" choice assures us that the spacetime has 
the usual black hole causal structure for v —^ ±oo and, 
consequently, that QNM can be defined in the usual way 
and that the corresponding frequencies can be properly 
compared. For all the cases considered, the inertial tran- 
sients dissipate away and (^(w) tends to follow the track 
of m~^{v) rather quickly, rendering the robustness of the 
numerical QNM analysis. The integration in double-null 
coordinates has turned out to be much more efficient than 
the integration in radiation coordinates'^, allowing us to 
reach the precision necessary to unveil the reported non- 
stationary behavior with quite modest computational re- 
sources. Despite the radiation coordinates are known to 
be defective at the horizon 0|, we believe that numerical 
analysis as those one presented in are still confident 
since the QNM analysis is always concerned to the exte- 
rior region of the black hole. However, as confirmed by 
the present calculations, algorithms based on the charac- 



teristic integration are typically by far more efficient. 

An interesting extension of this work would be analy- 
sis of the highly damped QNM (overtones, n > 0). Since 
for such overtones the ratio Iwi/wrI is always larger than 
for the n = QNM considered here, including, for suf- 
ficient large n, cases for which ja^i/wnl > 1, it would be 
interesting to check if the non-stationary inertial behav- 
ior could somehow be attenuated for n > 0. We notice 
that the numerical analysis presented here cannot be ex- 
tended directly to the n > case since one cannot iden- 
tify the n > frequencies with sufficient accuracy. We 
believe this could be attained, in principle, by means of 
the WKB approximation. 

Although the typical astrophysical situations of mass 
accretion for black holes will hardly keep spherical sym- 
metry intact during intermediate stages, our results can 
be used as a first approximation for the scattering by 
this sources since, as it is well known, after the tran- 
sient phases, the system should accommodate itself in a 
stationary spherical symmetric configuration. One must, 
however, keep in mind that non-stationary inertial be- 
havior in the QNM frequencies is expected to take place 
whenever \m"\ ^ jcjij, rendering the analysis of the rapid 
variating situations a subtler task. 

The evaporation by Hawking radiation could indeed 
be considered as a physically genuine process where the 
black hole mass decreases and spherical symmetry is 
maintained. Our approach can be applied for this case, 
with some crucial remarks about the causal structure of 
the spacetime left behind after the evaporation||23]. 

We finish by noticing, as it is well known, that the 
damped oscillations correspond to an intermediate phase 
in the wave scattering by asymptotic fiat black holes. 
The very last final phase corresponds, indeed, to a power 
law decay. In the problems considered here, the power 
law tail typically appears for large values of v for which 
the QNM already settled down in the stationary phase, 
with no trace from the transients. We do not detect, 
in this case, any influence of the time-dependence of the 
potential in these tails. There is no contradiction with 
the results reported in 0. This is a consequence of our 
choice of "asymptotic Schwarzschild" mass functions, for 
which the corresponding potential has no resemblance 
with those ones considered in 
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